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Abstract. We present new accurate CCD uvby light curves for the LMC eclipsing binaries HV982 and HV12578, 
and for the SMC systems HV1433 and HV11284 obtained at the Danish 1.5m telescope at ESO, La Silla. The 
light curves were derived from DoPHOT photometry, and typical accuracies are between 0.007 and 0.012 mag per 
point. Standard uvby indices have also been established for each binary, primarily for determination of interstellar 
reddening and absorption. For HV982 and HV12578, accurate photometric elements have been established. Both 
systems consist of two detached components of comparable sizes in an eccentric orbit. Adopting the spectroscopic 
elements given by Fitzpatrick et al. (J5002) for HV982, we derive absolute dimensions of its components which 
agree well with their results. A distance modulus of Vb — My = 18.63 ± 0.08 is obtained, corresponding to a 
distance of 52.6 ± 2.0 kpc, which is in formal agreement with (although slightly larger than) their determination. 
HV1433 and HV11284 both consist of two rather close, deformed and quite different stars. As the mass ratios 
between the components (and their rotation rates) are not known, definitive photometric elements can not yet 
be obtained, but we present a sample of possible photometric solutions. In a series of forthcoming papers we will 
combine our uvby observations with high-dispersion spectra from the UVES spectrograph on the ESO Very Large 
Telescope (VLT) and present absolute dimensions, chemical abundances and distances for selected LMC and SMC 
systems, including HV12578 and refined results for HV982. 
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binaries: close - Stars: binaries: eclipsing - Techniques: photometric 



1. Introduction 

Double-lined eclipsing binaries provide, apart from ob- 
jects where purely geometrical methods can be applied, 
the most direct distance determinations available. From 

Send offprint requests to: J.V. Clausen 

* Based on observations carried out at the Danish 1.5m tele- 
scope at ESO, La Silla, Chile 

** Tables 25-40 will be made available on electronic form 
at the CDS via anonymous ftp to 130.79.128.5 or via 
: / / cdsweb.u-strasbg.fr / Abstract .html 



complete, accurate observations and careful analyses of 
well-detached systems, distance moduli better than ±0.10 
mag can be obtained. 

With present-day facilities, including ground based 10- 
m class telescopes and the Hubble Space Telescope (HST) , 
precise studies of systems in Local Group galaxies, pri- 
marily the Magellanic Clouds (MC), have become possi- 
ble. Therefore, accurate LMC and SMC distances can now 
be established from eclipsing binaries, leading also to the 
calibrations needed for other important extragalactic di- 
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stance indicators, like e.g. Cepheids, reaching far beyond 
the Local Group. Reviews on the subject have recently 
been presented by Paczynski l|1997[) and Clausen ( 2000). 

Also, the accurate absolute dimensions obtained for 
the MC binary components makes it possible to gain in- 
sight into e.g. core structure (e.g Ribas et al I2000|l and 
mass loss for massive stars through comparison with stel- 
lar models. Furthermore, they provide fundamental mass- 
to- luminosity (M/L) relations for metal-deficient stars of 
e.g. typical LMC/SMC compositions. 

The first accurate CCD light curves of LMC and SMC 
eclipsing binaries were obtained at the Danish 1.5m tele- 
scope at ESO, La Silla by Jensen et al. I|1988(l . who give 
references to earlier photometry. The present paper should 
be seen as a continuation of this work. Meanwhile, de- 
dicated CCD light curve observations have been done for 
several LMC and SMC systems at Mt. John University 
Observatory, New Zealand (Pritchard et al. 11998a! I1998bl 
and references therein), and several thousands of new sy- 
stems have been discovered from the EROS, MACHO, 
and OGLE microlensing projects (e.g. Grison et al. 119951 
Alcock et al. 113971 Udalski et al. I1998bf) . 

Low to medium resolution spectroscopy is available 
for a few MC systems (e.g. Niemela & Bassino 1994), 
and moderately accurate distances and dimensions based 
on CCD light curves and radial velocities from such 
spectroscopy have been published for some systems, e.g. 
HV2241 (Pritchard et al. ll998aj) . HV2543 (Ostrov et al. 
IMKTjl . and HV5936 (Bell et al. IT§9l|) in the LMC, and 
HV1620 (Pritchard et al. Il998a| and HV2226 (Bell et al. 
1991) in the SMC. 

However, an important step forward was recently taken 
when a very accurate distance modulus (about ±0.07 
mag) and stellar dimensions (about ±2.5% for radii and 
±6% for masses) were determined for the LMC sys- 
tem HV2274 (Watson et al. IT5521 Guinan et al. IT5M1 
Udalski et al. 11998a! Nelson et al. 12113 Ribas et al. 123001 
Groenewegen & Salaris l2001[) . Masses and radii were ob- 
tained from normal analyses of light curves and radial ve- 
locity curves, based on rather few high dispersion spec- 
tra, whereas effective temperatures, reddenings and dis- 
tances were obtained from UV/optical spectrophotometry. 
Slightly deviating distances are obtained from the different 
studies, mainly due to differences in the reddening deter- 
minations, but they all tend to support the short LMC 
distance scale (e.g. Gibson [2000). 

Very recently, accurate distances and dimensions for 
two additional LMC systems, HV982 (Pritchard et al. 
I1998bl Fitzpatrick et al HUU^i and EROS1044 (Ribas et al. 
2002}, have become available. A comparison of the results 
for HV2274, HV982, and EROS1044, and a discussion of 
the derived distances to the optical center of the LMC is 
given by Ribas et al. (|2002|) . Clearly, studies of many more 
binaries, well distributed across the LMC, are needed, and 
are important also for investigations on the orientation 
and structure of the galaxy. SMC systems should also be 
included to extend the investigation to significantly lower 



Table 1. The observed eclipsing binaries. For comparison, 
the position of the optical centers of LMC and SMC are 
included. 



ID 




Pdays 


a(2000.0) 


5(2000.0) 


LMC 






05 ft 24 m 


-69°45' 


HV982 


14.7 


5.34 


05^29™ 53 s 


-69°09'23" 


HV12578 


15.8 


2.48 


05' l 21 m 32 s 


-66°21'15" 


SMC 






00^53™ 


-72° 50' 


HV1433 


16.4 


2.05 


00 h 47 m ll s 


-73°41'18" 


HV11284 


17.1 


3.63 


00 h 49 TO 43 s 


-72°51'10" 



metallicities, which will also allow a study of the metallic- 
ity dependence of the Cepheid period-luminosity relation. 

In this paper, we present new complete uvby light 
curves and analyses of them for the LMC systems HV982 
and HV12578, and the SMC system HV11284, and par- 
tial results for the SMC system HV1433. Studies of redde- 
ning and metallicity in the observed fields (about 6.5 x 6.5 
arcminutes) centered on the four systems have previously 
been published (Larsen 11996! Larsen et al. [2001 ) . 

Further, together with a group of colleagues we are 
undertaking a large scale project, combining light curve 
observations, standard photometry, and high dispersion 
spectroscopy from the ESO VLT UT2 telescope and the 
new very efficient UVES echelle spectrograph (Dekker et 
al. D'OdoricoEOOOl)- Until now, complete spectro- 

scopic observations have been obtained for four LMC sys- 
tems, including HV982 and HV12578, and two SMC sys- 
tems, and analyses are in progress. 

2. Selection of candidates 

When the uvby observations were planned, the thousands 
of newly discovered MC binaries from the microlensing 
projects were not yet available, and our source for the se- 
lection was then limited to the about 110 systems identi- 
fied mainly from the photographic material of the Harvard 
surveys (Shapley & Nail 1 19421 HUSH Payne-Gaposchkin & 
GaposhkinlTkHl G aposhk in IT5?ni Hodge & Wright 

ITT^TI [TT777I Kreiner[M2). 

About 10 candidates per galaxy were selected, but for 
several reasons (field crowding, faintness, difficult orbital 
periods, lack of observing time, etc.) the project was gra- 
dually concentrated on the four systems listed in Tabled 
Part of the fields, centered on the eclipsing binaries, are 
shown in Figs. □[213 andU 

HV982 is located about one degree NE of the opti- 
cal center of LMC, whereas HV12578 is located in the 
northern outskirts of the LMC Bar, about three degrees 
from the optical center. HV1433 is located in the sou- 
thern part of the central region of SMC, just south of 
the OGLE - II field SMC.SC4 (Udalski et al. H998b|l . and 
HV11284 slightly west of the optical center of SMC in the 
OGLE -II field SMC.SC5. HV11284 is identical to the 
OGLE - II eclipsing binary SMC.SC5.140701. 
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Fig. 1. HV982 and the surrounding field (1.5 arcmin 
square). The y image was observed at 1.0 arcsec seeing. 
North is up, and east is to the left. 




Fig. 2. HV12578 and the surrounding field (1.5 arcmin 
square). The b image was observed at 1.1 arcsec seeing. 
North is up, and east is to the left. 
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Fig. 3. HV1433 and the surrounding field (1.5 arcmin 
square). The y image was observed at 1.1 arcsec seeing. 
North is up, and east is to the left. 




Fig. 4. HV11284 and the surrounding field (1.5 arcmin 
square). The b image was observed at 1.1 arcsec seeing. 
North is up, and east is to the left. 



3. Observations and basic reductions 

The CCD observations were carried out during several pe- 
riods between November 1992 and November 1995 at the 
Danish 1.54 meter telescope at ESO, La Silla, equipped 
with a direct camera and CCD #28 (a thinned 1024 x 1024 
Tek device). This combination of telescope and CCD 
yields a field size of about 6.5 x 6.5 arcminutes. The im- 



age scale is 0.38"/pixel. The same set of uvby interference 
filters were used throughout, except for the early 1992 
b observations. Photometry from several b images of the 
same field obtained with the two filters were carefully com- 
pared, but no significant differences between the instru- 
mental systems were found. 

The binaries were always positioned near the field cen- 
ters, but no effort was made to place them on exactly 
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the same pixels on all exposures. Also, the orientation of 
the camera/chip varied between the different observing 
runs. Typical exposure times for HV12578, HV11284, and 
HV1433 were 5 (vby) and 20 minutes (u), but for the 
brighter HV982 we used 2.5 and 10 minutes, respectively. 
The observations were in general done in ybvyu sequences 
in order to obtain complete light curves in one band (y) 
and acceptable coverage in the other bands. 

Initial reductions, essentially dark current and bias 
subtractions and flatfielding, where carried out using stan- 
dard IRAF 1 tools. As the dark current was very low, only 
around 4 counts per hour with no 'hot spots' or other de- 
fects visible, we decided to use just a single common value 
for the dark current in order not to introduce additional 
noise. In each observing run, typically 10 bias exposures 
were made, and in all cases the bias level was found to 
be constant within 1-2 ADU during a given run, whereas 
the level changed slightly from year to year, depending on 
the electronic setup of the CCD. The bias frames were all 
very uniform, the only defect being a weak pattern from 
the thinning of the chip, and we used combined frames for 
the bias subtraction. In each observing run, sky flatfield 
exposures were made every evening and morning, and the 
individual frames, typically 5 per filter with slight offsets 
between them, were carefully checked against a reference 
frame obtained with the same filter. Any exposures that 
contained bright stars or other obvious anomalies were 
rejected. For the remaining exposures, we found no obvi- 
ous variations from night to night or between morning and 
evening flats, and we thus decided to use a single combined 
flatfield in each filter for each run. The accepted flatfield 
exposures were combined applying a sigma-clipping rejec- 
tion algorithm to avoid any (faint) stars that might be 
present in the images. 

As a last step, we measured the seeing from two or 
more stars in each frame, and determined the sky back- 
ground level. This information is given in Tables 25-40 
(light curve tables, instrumental system). No images with 
seeing above about 3 arcsec were accepted, and most im- 
ages were obtained at 1-1.5 arcsec seeing. Very few images 
with a background above 100 ADU were accepted. 

Additional information on the reduction is given by 
Larsen RTM5\i . 

4. Photometry and formation of light curves 

4.1. Light curves 

Several packages are available for crowded field pho- 
tometry, including MOMF (Kjeldsen & Frandsen 
Hnn2l, ROMAPHOT (Buonnano et al. IrMSjl . DoPHOT 
(Schechter et al. lT^I . and DAOPHOT ( Stetson IT§8T|) . 
Especially the latter two have been widely used, and while 
DAOPHOT may be the most flexible of the two, DoPHOT 

1 IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc. under contract 
with the National Science Foundation. 



Table 2. Mean DoPHOT errors (mag.). 



Binary 


y 


b 


V 


u 


HV982 


0.008 


0.007 


0.008 


0.008 


HV12578 


0.009 


0.008 


0.009 


0.010 


HV1433 


0.009 


0.008 


0.010 


0.012 


HV11284 


0.008 


0.007 


0.008 


0.009 



is particularly well suited for reductions of large quan- 
tities of data because only a minimum of user interac- 
tion is required. Tests of DoPHOT and DAOPHOT by 
ourselves and others (e.g. Larsen 1996; Ferrarese et al. 
2000) show that the two packages generally provide pho- 
tometry of comparable accuracy, as well as very reali- 
stic error estimates. Because of its large number of user- 
definable parameters and greater flexibility, DAOPHOT 
may give slightly more accurate photometry in difficult 
cases such as strongly crowded fields, or if the PSF varies 
significantly across the frame. This may be partly due 
to the more sophisticated PSF modeling in DAOPHOT. 
However, DoPHOT also produces very good photometry 
and runs substantially faster than DAOPHOT. Because 
of the large volume of data to be reduced in this project, 
we therefore opted to use DoPHOT for the analysis of the 
LMC and SMC fields. 

Following the DoPHOT analysis, the photometry out- 
put files were combined into a single database for each 
binary field using the DBMAG package developed by one 
of us (JS). Besides efficient book keeping and photometry 
extraction tools, it contains tasks for identification and de- 
tection of variable stars based on the technique by Welch 
& Stetson ([£993). For each binary field and photometric 
band, a master frame was selected, and all similar frames 
were then brought to the zero point of the master. 

As a final check of the procedure, we studied one of 
the databases, the y band of the HV11284 field, in further 
detail. A number of isolated, supposedly constant stars 
spanning the magnitude range of HV11284 through its 
orbit were extracted, and their mean magnitude and the 
corresponding r.m.s. per observation were calculated. No 
systematic patterns in the individual deviations from the 
means were seen for the stars, and r.m.s. values between 
0.009 and 0.011 mag were obtained. For all stars, the r.m.s. 
values are close to the DoPHOT error estimates of the 
single magnitudes, and we conclude that the frames have 
been safely brought to the same zero point. Therefore, we 
have chosen to apply directly the database results for the 
binary light curves, rather than basing them on a smaller 
number of preselected comparison stars, as is often done. 

The light curve observations in the instrumental sy- 
stem are presented in Tables 25-40, which will only be 
available in electronic form. For each light curve point, the 
internal frame number, the Heliocentric Julian Date corre- 
sponding to the time of mid-exposure, the orbital phase, 
the magnitude and its error (DoPHOT), the sky back- 
ground level (in ADU), and the seeing (in pixels, the scale 



J.V. Clausen et al.: uvby CCD light curves 



5 



is 0.38"/pixel) are given. Typical errors per point are close 
to 0.01 mag (Table 0). 

Phases were calculated from the individual 
ephemerides given below. When possible, times of 
minima were calculated from the Kwee & Van Woerden 
( 1956 ) method. Ephemerides were, if possible, determined 
from least squares fits to new and published times of 
minima, else by fitting the observations together applying 
the method by Lafler & Kinman (1965). The resulting 
light curves are shown in Figs. El El El and 1161 



5. Standard uvby indices 

In order to obtain reliable apparent magnitudes and co- 
lors for the eclipsing binaries, essential for reddening and 
distance determinations, a careful transformation to the 
uvby standard system is needed. We have based this cali- 
bration on numerous CCD observations of secondary uvby 
standard, for which photoelectric observations were simul- 
taneously done at the Stromgren Automatic Telescope 
(SAT) at ESO, La Silla. The SAT observations also pro- 
vided accurate nightly extinction coefficient needed for the 
calibration of the CCD photometry. The SAT observations 
and standard uvby indices for the secondary standards are 
described by Clausen et al. 1)1997(1 . 

Because DBMAG automatically adjusts the photom- 
etry of all individual frames to a common zero-point, it 
was sufficient to derive a photometric transformation to 
the uvby standard system for the reference frames. For 
this purpose, we first used the substar task in the IRAF 
DAOPHOT package to subtract all stars except a few 
bright ones from the reference frames. Aperture photome- 
try for these stars was then obtained with the phot task, 
thus providing a set of "tertiary standard stars" in each 
frame (see Larsen et al. 2001] for details). Finally, the 
DoPHOT photometry was tied into the standard system 
by determining the offsets between the DoPHOT magni- 
tudes and aperture photometry for the tertiary standard 
stars. Figure0shows the difference between the DoPHOT 
and aperture photometry for the selected tertiary stan- 
dards in the HV 982 field. Typical uncertainties on the 
zero-points were estimated to be 0.01 - 0.02 mag. 

We further checked the accuracy of the photometric 
zero-points by comparing DAOPHOT PSF-fitting pho- 
tometry on the co-added frames (Larsen et al. I2()01|l 
with the mean DoPHOT magnitudes from the DBMAG 
database. For the LMC fields, the comparison showed ex- 
cellent agreement between the DoPHOT and DAOPHOT 
photometry, with systematic differences < 0.01 mag in y 
and b — y and < 0.02 mag in mi and c\. For the SMC 
fields we found slightly larger differences of ~ 0.02 mag in 
y and up to 0.03-0.04 mag in mi, but still compatible with 
the estimated zero-point uncertainties on the photometry. 
Note that the zero-points were established for the indi- 
vidual passbands, leading to larger formal uncertainties 
on indices such as mi and ci, which are based on three 
bands each. 
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Fig. 5. Zero-point differences between aperture and 
DoPHOT photometry for "tertiary standards" in the 
HV982 field. The dashed lines indicate the adopted zero- 
points. 



The resulting standard uvby indices for the eclipsing 
binaries are given in Table [21 Intrinsic standard uvby in- 
dices, interstellar reddening, and absorption derived from 
the (b — y)o — Co relation by Crawford l|1978[) are given in 
Table 0] For HV982, the reddening determined from the 
uvby photometry agrees well with that recently published 
by Fitzpatrick et al. lf2UU2|l . E(B - V) = 0.086 ± 0.005, 
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Table 3. Standard uvby indices for the eclipsing binaries (no correction for reddening applied) at maximum light level 
outside eclipses. N is the number of observations used to form the mean value, a is the r.m.s. error (per mag/index) 
given in mmag. 



ID 


V 


a 


N 


(b-y) 


a 


N 


mi 


a 


N 


Cl 


a 


N 


HV982 


14.648 


8 


69 


-0.042 


9 


47 


0.070 


15 


44 


0.020 


17 


42 


HV12578 


15.768 


6 


27 


-0.050 


10 


20 


0.094 


13 


16 


0.077 


25 


11 


HV1433 


15.547 


11 


3 


-0.070 


5 


4 


0.106 


14 


4 


0.014 


25 


3 


HV11284 


15.179 


7 


15 


-0.023 


10 


27 


0.094 


15 


21 


0.069 


13 


17 



Table 4. Intrinsic standard uvby indices at maximum light level outside eclipses, and interstellar reddening and 
absorption for the eclipsing binaries. 



ID 


V 


(b-y) 


m 


CO 


E(b-y) 


E(B - V) 


A v 


HV982 


14.350 


-0.111 


0.091 


0.006 


0.069 


0.097 


0.298 


HV12578 


15.529 


-0.106 


0.111 


0.066 


0.056 


0.078 


0.239 


HV1433 


15.369 


-0.111 


0.118 


0.006 


0.041 


0.057 


0.178 


HV11284 


14.818 


-0.107 


0.119 


0.052 


0.084 


0.118 


0.361 



based on a completely different and much more complex 
approach. 

6. HV982 (LMC) 

Already the photographic Harvard Survey data 
(Gaposhkin I1970|l show that HV982 is a well-detached 
eclipsing binary, presumably with an eccentric orbit. CCD 
light curves obtained at Mt. John University Observatory, 
New Zealand by Pritchard et al. fl998bl clearly confirm 
these characteristics. Furthermore, with new minima 
observations, apsidal motion with a period of about 200 
yr could be detected, making HV982 just the second MC 
eclipsing binary with measured orbital motion (the first 
known case is HV2274, also in the LMC; Watson et al. 
ITMSl Ribas et al. I2UUU|). 

Pritchard et al. (1998b) also did standard uvby pho- 
tometry and IUE ultraviolet spectroscopy, and deter- 
mined interstellar reddening and photometric elements for 
HV982. Subsequently, Pritchard & Niemela l|2UUU|l pre- 
sented radial velocities and derived mean masses and radii 
for the components, and a distance modulus of Vq — My = 
18.7 ± 0.15 was obtained. Very recently, Fitzpatrick et 
al. I|2002(l published accurate dimensions for HV982 and 
found its distance modulus to be Vo — My — 18.50 ± 0.06. 
The results were obtained by adding HST UV/optical 
spectrophotometry and radial velocities based on high- 
resolution optical spectroscopy from the Blanco 4-m tele- 
scope at Cerro Tololo, Chile to existing photometry for 
HV982. 

As HV982 is clearly an important LMC eclipsing bi- 
nary, we decided - parallel to the Mt. John observations 
- to gain further light curve data. Mainly because the 
field near the system is quite complicated, with several 
close neighbours (Fig. P) , a delicate subtraction procedure 
was needed at the average Mt. John seeing of 3-5 arcsec 
(Pritchard et al. 1998b). This can be avoided at the ave- 



rage La Silla seeing of about 1-2 arcsec. Also, better uvby 
standard indices were needed. 

In this paper, we present the light curves and photo- 
metric elements derived from them, plus standard uvby 
photometry. Also, absolute dimensions and distances are 
derived, but we expect to publish significantly improved 
results when the analyses of the VLT/UVES spectra has 
been completed. 

6.1. Light curves 

HV982 was observed on 43 nights between JD2448970 and 
JD2450047, and 314(y),124 (b), 115 (v), and 103 (u) points 
were obtained in the four bands. The light curves are given 
in Tables 25-28 and shown in Fig. |SJ The coverage in y, 
b, and v is fairly complete, whereas a gap exists in the 
u observations of primary eclipse. Typical accuracies per 
point are given in Table [3 

The two eclipses are of comparable depths with central 
secondary eclipse occurring at phase 0.428. No significant 
change of the (b—y) color is seen during the eclipses. Slight 
reflection effects are noticed outside the eclipses. 

6.2. The close companions 

As seen on Fig. ^ HV982 has several rather close compa- 
nions. When possible, Pritchard et al. I|1998b|l PSF fitted 
as a group HV982, the two brightest companions at about 
3 arcsec distance (of which the southern actually seems 
to be two very close stars), and a third at about 6 arcsec 
distance, except at too large seeing where the contribu- 
tions of the companions were subtracted assuming known 
positions and magnitudes. 

As our observations were done at much better seeing, 
we have been able to isolate HV982 and avoid signifi- 
cant contamination from all the companions in the PSF 
fits to HV982. This is demonstrated in Fig. Here the 
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Fig. 6. uvby light curves of HV982 (LMC). The theoretical curves correspond to the solutions given in Table [5] 



database photometry for the 3 arcsec distance companion 
(id2673) NNE to HV982 is plotted versus the orbital phase 
of HV982 at the time of the observation, versus the see- 
ing, and versus the sky background level. The companion, 
which is the brightest one, is about 3.4 mag (y) fainter 
than HV982 outside eclipses and its DoPHOT photom- 
etry therefore much less accurate. However, we find no 
correlation between the magnitude of the companion and 
the orbital phase (i.e. light level) of HV982 or the sky 
background level. There is, perhaps, a tendency of slightly 
fainter companion magnitudes for seeing above about 1.6 
arcsec where the scatter on the other hand also increases 
significantly. 

6.3. Ephemeris and apsidal motion 

One time of each of primary and secondary minimum, 
derived from the uvby CCD observations, are given in 
Table |S] together with times published by Gaposhkin 
lfT§7?|> and Pritchard et al. l|1998bj) . 

Apsidal motion parameters were derived from a 
weighted least squares method, following the formalism by 
Gimenez & Garcia-Palayo fTj5H3} and Gimenez & Bastero 




-6000 -4000 -2000 
Epoch 

Fig. 8. Apsidal motion for HV982 (LMC). O-C are the 
residuals (days) from the linear part of the apsidal mo- 
tion ephemeris defined by the parameters given in Table 
The full and dashed curves represent predictions from 
the apsidal motion parameters for primary (triangles) and 
secondary (squares) eclipses, respectively. 



(1995J). The orbital inclination i and eccentricity e were 
fixed to the values derived from the photometric analysis; 
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-0.2 0.0 0.2 0.4 0.6 0.8 -0.2 0.0 0.2 0.4 0.6 0.S -0.2 0.0 0.2 0.4 0.6 0.8 -0.2 0.0 0.2 0.4 0.6 0.8 

Phase Phase Phase Phase 



2673 y 2673 b 2673 v 2673 u 




100 200 300 400 100 200 300 400 100 200 300 400 100 200 300 400 

Sky background (ADU) Sky background (ADU) Sky background (ADU) Sky background (ADU) 



Fig. 7. DoPHOT photometry for the close companion (id2673) about 3 arcsec NNE to HV982 (LMC); see Fig. □ 
id2673 is about 3.4 mag fainter than HV982. The full drawn horizontal lines represent the mean magnitudes, and 
dashed lines are drawn at the 1 sigma levels. 



see Sect. 16.41 The results, which are shown in Tabic Eland 
Fig. |S1 agree quite well with those obtained by Pritchard 
et al. I|1998bjl . 

For the calculation of phases of the uvby light curves, 
we adopt the following linear ephemeris derived for the 
mid-epoch of the observations from the apsidal motion 
parameters: 



Min I = 2449340.7172 + 5 d 335220 x E (1) 



The adopted period is identical to the sidereal period. 
During the light curve observations from JD2448970 to 
JD2450047, the longitude of periastron has changed from 
224?5 to 225?0 and the phase of mid secondary eclipse 
relative to mid primary eclipse from 0.428 to 0.429. For 
the light curve analyses, the effect of apsidal motion can be 
ignored, except for a few points during secondary eclipse 
from JD2450041, which have been excluded. 



6.4. Light curve analysis, photometric elements 

Pritchard et al. l)1998bj) presented photometric elements 
derived from their CCD light curves, adopting a modi- 
fied version of the Wilson-Devinney (WD) code (Wilson 
& Devinnev ll971l Wilson[l994). Light curves in all bands 
were analyzed simultaneously, and four possible sets of 
elements were presented, as the differential least squares 
fit did not converge towards one unique solution. This re- 
flects a situation which is often seen for partially eclipsing 
systems with fairly identical components, both for circular 
and eccentric orbits. It is due to significant correlation be- 
tween several elements. Often the difference between the- 
oretical light curves fitted to the observations for a large 
range of adopted k — r s /r p (p referring to the compo- 
nent eclipsed at the deeper primary eclipse at phase 0.0) 
are identical within a few mmag. Pritchard et al. 0"998b) 
furthermore did not know the mass ratio q between the 
components, which is needed for calculation of the stellar 
shapes. The classical ways to break (part of) the degener- 
acy are to include independent spectroscopic information 
on the luminosity ratio (see e.g. Andersen et al. 11983)1 . or 
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Table 5. Times of minima for HV982 (LMC). 0-C values 
(in days) are calculated for the apsidal motion parame- 
ters given in Table References are: G1977, Gaposhkin 
IWTI P1998, Pritchard et al. I1998hl C2003, this pa- 
per. The r.m.s. errors of the G1977 times are estimated. 
Note the perfect agreement between P1998 and C2003 for 
JD2449337. 



HJD-2400000 


r.m.s. 


Type 


O-C 


Reference 


13946.555 


0.050 


P 


0.0290 


G1977 


17590.584 


0.050 


P 


0.0350 




23875.527 


0.050 


P 


-0.0360 




25849.645 


0.050 


P 


0.0130 




26060.243 


0.100 


s 


0.0900 




26412.253 


0.050 


s 


-0.0160 




26577.631 


0.050 


s 


-0.0270 




27786.315 


0.050 


p 


-0.0380 




29189.469 


0.100 


p 


-0.0700 




29629.338 


0.050 


s 


-0.0060 




31304.630 


0.050 


s 


0.0520 




32070.603 


0.050 


p 


0.0002 




33153.625 


0.050 


p 


-0.0420 




49335.3866 


0.0004 


p 


0.0018 


P1998 


49337.6668 


0.0004 


s 


-0.0002 




49337.6670 


0.0010 


s 


0.0018 


C2003 


49340.7172 


0.0005 


p 


-0.0028 




50695.8520 


0.0110 


p 


0.0086 


P1998 



Table 6. Apsidal motion parameters for HV982 (LMC). 



Parameter 


Value and r.m.s. error 


i n 


88.7 (assumed) 


e 


0.159 (assumed) 


To 


2449335.17575 ± 0.00049 


P anomalistic (d) 


5.335595 ± 0.000025 


Psidereal (d) 


5.335220 


wo n 


224.67 ±0.15 


lui (° / cycle) 


0.00253 ± 0.00018 


U(yr) ' 


208 ± 15 



- for eccentric systems - independent information on the 
orbital eccentricity from apsidal motion analyses (see e.g. 
Clausen et al. H^jl . 

Fitzpatrick et al. H2002[l determined q = M s /M p = 
1.029 ± 0.027 and the spectroscopic luminosity ratio 
L s /Lp = 1.15±0.05 (blue) between the components which, 
although the accuracy is not as high as could be desired, 
enabled them to exclude three of the solutions presented 
by Pritchard et al. <|1998bjl . 

In the following, we will present photometric elements 
derived from WINK (Wood 119711 extensions from Vaz 

rnreirresfii v az & Nordinnd rrwn N ordiun d & va z mm 

and WD (Wilson & Devinnev 113711 Wilson IT^Mjl . exten- 
sions from Vaz et al. I1995fl analyses of the new uvby light 
curves. Both codes provide physical binary models, ade- 
quate for the description of HV982. They differ mainly 
with respect to the model of the stellar shapes (3-axial 



Table 7. WINK solutions (y) for HV982 for a range of 
fixed k values between 0.90 and 1.10. T p = 24200 K and 
x = 0.33 have been assumed. 



y v y y v 





88.71 


88.26 


88.15 


88.25 


88.66 


±14 


±10 


±10 


±11 


±16 


ecostj 


-0.1115 


-0.1115 


-0.1115 


-0.1116 


-0.1117 




±2 


±2 


±2 


±2 


±2 


esintj 


-0.1243 


-0.1194 


-0.1170 


-0.1137 


-0.1124 




±43 


±43 


±43 


±44 


±44 




0.1669 


0.1634 


0.1617 


0.1593 


0.1585 


In 


228.11 


226.96 


226.38 


225.53 


225.20 


k 


0.90 


0.95 


1.00 


1.05 


1.10 


r p 


0.2138 


0.2085 


0.2033 


0.1982 


0.1931 




±7 


±7 


±7 


±7 


±7 


r a 


0.1924 


0.1981 


0.2033 


0.2081 


0.2125 


r v + r s 


0.4062 


0.4066 


0.4066 


0.4063 


0.4056 


Ts (K) 


23740 


23613 


23572 


23577 


23650 


±103 


±104 


±101 


±104 


±104 


J s 1 Jp 


0.959 


0.950 


0.947 


0.947 


0.952 


L s /Lp 


0.773 


0.856 


0.948 


1.049 


1.160 


cr(mag) 


0.0098 


0.0099 


0.0098 


0.0100 


0.0100 



Table 8. WINK solutions for HV982. k = r s /r p = 1.10 
and T v — 24200 K has been assumed. 





y 


b 


V 


u 


in 


88.66 


88.96 


88.84 


88.80 




±16 


±32 


±26 


±28 


ecosu 


-0.1117 


-0.1106 


-0.1106 


-0.1103 




±2 


±4 


±4 


±5 


esino; 


-0.1124 


-0.1102 


-0.1065 


-0.1097 




±44 


±75 


±71 


±76 


e 


0.1585 


0.1561 


0.1535 


0.1555 


<*n 


225.20 


224.90 


223.91 


224.83 


r p 


0.1931 


0.1929 


0.1912 


0.1932 




±7 


±12 


±11 


±12 


r s 


0.2125 


0.2122 


0.2103 


0.2125 


r p + r s 


0.4056 


0.4051 


0.4015 


0.4057 


Xp — x s 


0.33 


0.36 


0.39 


0.39 


Ts (K) 


23650 


23671 


23420 


23468 


±104 


±150 


±138 


±107 


J s I Jp 


0.953 


0.954 


0.931 


0.907 


L s 1 Lp 


1.160 


1.160 


1.134 


1.105 


cr (mag) 


0.0100 


0.0104 


0.0097 


0.0100 



ellipsoids with 4th and 5th order terms included versus 
Roche potentials) and with respect to the elements de- 
termined through the differential least squares procedure 
(e.g. r p and k versus potentials, and ecosw and esinw ver- 
sus e and lo). WINK is faster and more used-friendly than 
WD and therefore well-suited for detailed investigations 
of multi-dimensional parameter spaces. An extended IDL- 
FORTRAN WINK version has been developed by one of 
us (JVC) for that purpose. 

We have selected to perform independent rather than 
combined analyses of the uvby light curves in order to 
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obtain full understanding of the accuracy to which the 
elements can be established. 

6.4.1. WINK analyses 

The results from a series of WINK analyses for fixed k val- 
ues between 0.88 and 1.12 are shown in Fig. El and listed 
in Tabled An effective temperature of T p = 24200 K was 
assumed in agreement with the combined intrinsic uvby in- 
dices given in Tableland the calibration by Napiwotzki 
et al. (|1992f) . as well as with the energy distribution ana- 
lysis by Fitzpatrick et al. (2002). Surface fluxes for the 
uvby bands, as function of T c g and \ogg were based on 
Kurucz (1992) ATLAS9 atmosphere models (although for 
Solar composition). Linear limb darkening coefficients x 
(see TablelSJ) by Diaz-Cordoves et al. (|1995[) were adopted, 
and gravity darkening exponents of 0.25 and bolometric 
albedos of 1.0 were assumed in accordance with atmo- 
spheres in radiative equilibrium. For the model calcula- 
tion of the stellar deformations, a mass ratio of q = 1.0 
was assumed, and the rotation of the stars was assumed 
to be pseudosynchronized (rate at periastron). Initially it 
was settled that no phase shift was needed, i.e. the epoch 
of the ephemeris given in Eq. JIJ is well established. 

As seen from the upper panel of Fig. [51 solutions for k 
between 0.90 and 1.10 fit the observations nearly equally 
well in all four bands. Several parameters correlate directly 
with k, primarily e (through esimu) and L s /L p , whereas 
T s (and thereby the ratio of surface fluxes J s / J p ) and the 
sum of relative radii are nearly constant. Solutions from 
the four bands agree quite well - the y solution carries 
highest weight due to the larger number of points and the 
most complete phase coverage for this band. A luminosi- 
ty ratio (vby) of L s /L p — 1.15 ± 0.05, as determined by 
Fitzpatrick et al. lj2l)0"2Jl . corresponds to k = 1.10 ± 0.03. 
We have assumed this value, and solutions for k = 1.10 
are presented in Table 03 and (O — C) residuals (observed 
minus calculated mag) are shown in Fig. 1101 No systematic 
effects with phase and seeing are seen, but for an unknown 
reason a slight tendency of positive residuals (observation 
fainter than theoretical point) for background levels above 
about 150 ADU seems to be present. However, the r.m.s. 
errors of the fits are comparable to the photometric errors 
of the CCD magnitude differences. 

6.4.2. Wilson-Devinney analyses 

For comparison, both with the results from our WINK 
analyses and with the results by Fitzpatrick et al. ( 2002 ) , 
we present in Table H3 a series of WD solutions for the y 
light curve. The WD solutions were done for very nearly 
the same ratio of radii k as the WINK solutions given 
in Table [7| and as for the WINK solutions, a mass ratio 
of q = 1.0 and pseudosynchronized rotation of the com- 
ponents have been assumed, y surface fluxes were based 
on Kurucz Q1992[l ATLAS9 atmosphere models, a linear 
limb darkening law with x = 0.33, gravity exponents of 



Table 9. WD solutions (y) for HV982 for fixed Q p . 
T p = 24200AT and x = 0.33 have been assumed. f2 specify 
surface potentials and r relative volume radii. 





V 


y 


y 


y 


y 




88.71 


88.34 


88.26 


88.32 


88.74 


±9 


±7 


±5 


±6 


±11 


e 


0.1633 


0.1628 


0.1606 


0.1593 


0.1594 




±20 


±21 


±19 


±19 


±22 


w ( ) 


2,1b. lb 


220.08 


oo c on 
220. 80 


2,20.64: 


220.3/ 


±67 


±72 


±68 


±69 


±79 




5.900 


6.020 


6.125 


6.240 


6.385 


n s 


6.488 


6.362 


6.253 


6.108 


6.015 




±18 


±18 


±16 


±17 


±16 


r p 


0.2172 


0.2115 


0.2066 


0.2016 


0.1958 


r a 


0.1958 


0.2007 


0.2061 


0.2112 


0.2155 


k 


0.902 


0.949 


0.998 


1.048 


1.101 


r p + r a 


0.4130 


0.4122 


0.4127 


0.4128 


0.4113 


Ts (K) 


23536 


23498 


23471 


23452 


23520 


±75 


±83 


±73 


±65 


±85 


L s /Lp 


0.778 


0.859 


0.947 


1.044 


1.157 


cr(mag) 


0.0097 


0.0097 


0.0097 


0.0098 


0.0099 



1.0 (which correspond to input values of 0.25 for WINK), 
and reflection albedos of 1.0, were adopted. 

When comparing the results from the two codes, it 
should first be noted that in WINK, the relative radii 
are corrected for the (small) expansion effects from ro- 
tation and gravitational interaction between the compo- 
nents, whereas such corrections are not applied in the WD 
code. In the case of HV982, the corrections are about 
0.002 for both components; see Table EH for further de- 
tails. Taking this code difference into account, the two 
binary models are seen to produce practically identical 
results for HV982. 

6.4.3. Adopted photometric elements 

From the results of the uvby light curve analyses de- 
scribed above, we derive the mean photometric elements 
for HV982 listed in Table [IJB Within the errors, our re- 
sults agree well with those given by Fitzpatrick et al. 2002 
Slightly more accurate values have been obtained for the 
orbital inclination and eccentricity, probably because our 
uvby photometry is more accurate than that by Pritchard 
et al. ©98b). Accepting the luminosity ratio derived by 
Fitzpatrick et al. ( 2002J) , both analyses yield relative radii 
accurate to about 1%, i.e. at the level of the expansion of 
the components compared to single stars; see Sect. 16.4.21 

6.5. Absolute dimensions, reddening, and distance 

Preliminary absolute dimensions and distance for HV982 
are given in Table ITT1 together with photometric data for 
its components. The results are based on the mean pho- 
tometric elements listed in Table ITUl and the spectroscopic 
elements by Fitzpatrick et al. 120021 The uvby indices have 
been derived from the combined indices at phase 0.25 
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90.0 
89.5 
89.0 
88.5 
88.0 
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0.70 
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k 
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0.90 0.95 1.00 1.05 

k 
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0.01 10 ! 
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0.0080 E 
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1 
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- 
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0.90 0.95 1.00 
k 




0.90 0.95 1.00 
k 




0.95 1.00 1.05 1.10 
k 




0.90 0.95 1.00 1.05 

k 




0.90 0.95 1.00 1.05 1.10 
k 



0.4070 

0.4060 | 

0.4050 

0.4040 

0.4030 | 

0.4020 

0.4010 




0.90 0.95 1.00 1.05 1.10 
k 




2.34x10 



0.90 O.i 



1. 00 1.05 1.10 
k 



Fig. 9. Results from WINK analyses of HV982 for adopted k = r s /r p values between 0.88 and 1.12. The individual 
figures show s: r.m.s. error (mag), si: r.m.s. error (mag) for primary eclipse, s2: r.m.s. error (mag) for secondary 
eclipse, i: orbital inclination i, e: orbital eccentricity e, w: longitude of periastron to (°), rl: relative radius for primary 
component r p , r2: relative radius for secondary component r s , rl2: sum of relative radii, 121: luminosity ratio L s /L p , 
jn: surface flux ratio, t2: effective temperature of secondary component T s (K). Symbols are: cross y, diamond &, 
triangle v, square u. 



(Table|3) and the luminosity ratios ( Table IT0|) . Reddening, 
interstellar absorption and intrinsic colors are based on 
the calibration by Crawford (1978). Effective tempera- 
tures were calculated from the empirical [u — b] calibration 
by Napiwotzki et al. Ijl992l) . and the well established tem- 
perature difference from the light curve solution, and bolo- 
metric corrections are based on the calibration by Flower 

The dimensions and temperatures we obtain for 
HV982 agree well with those presented by Fitzpatrick et 
al. |2)02 , although the photometric data and the methods 
used for the analyses are different. The distance has been 
derived from the 'classical' formula 

V - M v = V - A v + 5 x log(i?/i? )+ 

10 x log(T cff /T cff0 ) + BC- Afboio 1 ' 

adopting r.m.s. errors of 0.01 (V), 0.02 {Ay), and 0.05 
(B.C.). Our value is slightly larger, but in formal agree- 



ment with that obtained by Fitzpatrick et al. 120021 f50.2± 
1.2 kpc). However, we expect to refine both absolute di- 
mensions and distance when independent spectroscopic 
elements, temperatures, luminosity ratios, and chemical 
abundances become available from the VLT/UVES spec- 
tra mentioned in Sect-^ 

Finally, we want briefly to comment on an interes- 
ting empirical method presented recently by Salaris and 
Groenewegen l|2002fl . The method, which is based on a 
calibration of Stromgren photometry, can be used for 
estimation of distances to B-stars in eclipsing binaries, 
and among various systems, they applied the method to 
HV982. For both components, a distance of 48.3 ±5.4 kpc, 
close to the result of the binary analysis, was obtained. It is 
based on the photometry by Pritchard et al. I|1998a|l . and 
the radii and V data from an earlier version of Fitzpatrick 
et al. 120021 which was, however, later revised significantly. 
We have repeated the calculations using the new results 
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HV982 y 




0.00 
-0.02 




Fig. 10. O — C residuals for the WINK solutions of HV982 for k = 1.10 given in Table |S1 plotted versus orbital phase, 
seeing, sky background, and DoPHOT error (mag). 



listed in Table and derive distances of 56.5 kpc and 
70.9 kpc, respectively, for the components of HV982. Such 
strongly deviating results might indicate that the calibra- 
tion by Salaris and Groenewegen needs further refinement. 



7. HV12578 (LMC) 

The photographic light curves of HV12578 (Gaposhkin 
11971) fl indicate a detached system in a circular orbit with 
components of comparable surface fluxes. The system is 
relatively bright, m pg = 15.05 at maximum light, and 
both eclipses are about 0.7 mag deep. The orbital pe- 
riod is 2*?48, and from the widths of the eclipses, relative 
radii of the components of about 0.20-0.25 are estimated. 
Although reflection effects seem to be present, we judged 
from the photographic data that HV12578 could very well 
be a promising candidates for determination of accurate 
dimensions and distances. The field near the binary is not 
much crowded, but our first CCD exposures showed that 
HV12578 has two much fainter neighbours at about 3 arc- 
sec distances; see Fig. |2 At average La Silla seeing they 



do not cause any problems for the PSF photometry of 
HV12578. 

7.1. Light curves 

We observed HV12578 on 33 nights between JD2448970 
and JD2450040, and 194 (y), 71 (6), 65 (v), and 53 (u) 
points were obtained in the four bands. The light curves 
are given m Tables 29-32 and shown in Fig. El Although 
the phase coverage is not as complete as we would prefer, 
the y and b light curves are adequate for photometric ana- 
lyses. The v and u data are more scarce and, besides pro- 
viding luminosity ratios between the components in these 
bands, they mainly serve as check data. Typical accuracies 
per point are given in Table [5] 

The new CCD observations reveal that reflection ef- 
fects are smaller than indicated by the Harvard material, 
and also that secondary eclipse is shifted slightly from 
phase 0.50, i.e. the orbit is eccentric. Both eclipses are 
deep, about 0.8 mag in all four bands, so the orbital incli- 
nation must be very close to 90° (a similar galactic F-type 
case is DM Vir; see Andersen et al. H984[) . From the CCD 
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Fig. 11. uvby light curves of HV12578 (LMC). The theoretical curves correspond to the solutions (k — 1.00) given in 
Table U31 



observations alone, it is not possible to settle definitively 
which eclipse is the deeper one. We have therefore decided 
to keep the To definition given by the Harvard ephcmeris 
fPavne-Gaposchkin ll971)) . although the final photometric 
analysis presented in Sect. 17.31 shows that for this defini- 
tion, the component eclipsed near phase 0.49 is marginally 
hotter than the other. We will return to this matter once 
the mass ratio becomes available from the analysis of the 
VLT/UVES spectra. 



7.2. Ephemeris 

Three times of primary minimum (at phase 0.00) and one 
of secondary minimum (near phase 0.49), derived from the 
uvby CCD observations, are given in Table El together 
with times published by Gaposhkin l|1970fl . Weighted 
(weight proportional to square of inverse r.m.s.) least 
squares fits to the times of primary and secondary minima 
yield periods of 2<?4803769 ± 0.0000017 and 2^4803730 ± 
0.0000017, respectively, and nearly identical results are 



obtained if equal weights are applied. For the calculation 
of phases of the uvby observations, we adopt 

Min 1 = 2448971.7778 + 2 d 4803769 x E , , 
±6 ±17 {A) 

The orbital period agrees well with that published by 
Payne-Gaposchkin 1)1971(1 . 2.48038, but is much more ac- 
curate. The time of secondary eclipse determined from the 
uvby observations reveals that the orbit of HV12578 is 
slightly eccentric, and this is confirmed by the light curve 
analysis; see Sect. 17.31 The slight difference in the peri- 
ods derived from times of primary and secondary eclipses, 
respectively, could in principle be due to apsidal motion, 
but since the orbital eccentricity is just about 0.015 this 
is not very likely, and it can not at all be confirmed from 
an analysis of the available times of minima. 

7.3. Light curve analysis, photometric elements 

For the analyses of the CCD light curves of HV12578, 
we have adopted the WINK and Wilson-Devinney models 
described in Sect. 16.41 
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Table 


13. WINK solutions for HV12578 for k = rj 


r p values of 0.98, 1.00, and 1.02. 


z = 90?0 and T p 


= 22000 


were adopted. For the b, v 


, and u 


solutions, 


ecostu and 


esincj 


were fixed at the y 


results. 










y 


b 


V 


u 


V 


6 


V 


u 


y 


b 


V 


u 


k 


0.98 


0.98 


0.98 


0.98 


1.00 


1.00 


1.00 


1.00 


1.02 


1.02 


1.02 


1.02 


esinu; 


0.00648 
±551 


0.00648 


0.00648 


0.00648 


0.00745 
±546 


0.00745 


0.00745 


0.00745 


0.01076 

±555 


0.01076 


0.01076 


U.UIUIO 


ecosw 


-0.01341 

±31 


-0.01341 


-0.01341 


-0.01341 


-0.01341 
±30 


0.01341 


-0.01341 


-0.01341 


-0.01341 
±31 


-0.01341 


-0.01341 


-0.01341 




0.0149 


0.0149 


0.0149 


0.0149 


0.0153 


0.0153 


0.0153 


0.0153 


0.0172 


0.0172 


0.0172 


0.0172 




154.2 


154.2 


154.2 


154.2 


150.9 


150.9 


150.9 


150.9 


141.3 


141.3 


141.3 


141.3 


Tp 


0.2450 
±10 


0.2479 
±13 


0.2442 

±2 


0.2497 
±25 


0.2426 
±10 


0.2455 

±14 


0.2419 

±20 


0.2472 
±23 


0.2402 
±10 


0.2430 

±13 


0.2395 

±22 


0.2447 
±23 


r B 


0.2401 


0.2429 


0.2394 


0.2447 


0.2426 


0.2455 


0.2419 


0.2472 


0.2450 


0.2478 


0.2443 


0.2496 


Tp + T s 


0.4851 


0.4908 


0.4836 


0.4944 


0.4852 


0.4910 


0.4838 


0.4944 


0.4852 


0.4908 


0.4838 


0.4943 




0.32 


0.37 


0.39 


0.39 


0.32 


0.37 


0.39 


0.39 


0.32 


0.37 


0.39 


0.39 


T s (K)' 


22328 
±76 


22328 
±100 


22447 
±151 


22079 
±107 


22332 
±75 


22295 
±100 


22430 
±151 


22073 
±108 


22277 
±76 


22289 
±100 


22409 
±151 


22058 
±108 


J a /J P 


1.027 


1.028 


1.041 


1.010 


1.027 


1.025 


1.039 


1.009 


1.023 


1.024 


1.037 


1.007 


L s /L p 


0.986 


0.987 


0.999 


0.971 


1.028 


1.026 


1.040 


1.011 


1.066 


1.068 


1.081 


1.051 


o-(mag) 


0.0109 


0.0093 


0.0148 


0.0157 


0.0107 


0.0094 


0.0149 


0.0159 


0.0111 


0.0093 


0.0149 


0.0160 



Table 10. Adopted photometric elements for HV982. r p 
and r s are corrected for expansion due to rotation and 
tidal effects, the ro results are not. The individual flux 
and luminosity ratios are based on the mean stellar and 
orbital parameters. 



Table 11. Astrophysical data for HV982. 
M bol Q — 4.74 has been assumed (Flower 1996). 



i 


88?7±0?2 


ecosiv 


-0.1117 ±0.0002 (at JD 2449500) 


esino; 


-0.1124 ±0.0044 (at JD 2449500) 


e 


0.159 ±0.003 


us 


225? 2 ± 1?1 (at JD2449500) 


r p 


0.194 ±0.003 


r„ 


0.214 ±0.003 


r p ± r s 


0.408 ± 0.003 


k 


1.10 ±0.03 (assumed) 




y b v u 


J s j Jp 


0.956 0.954 0.939 0.912 


L s j Lp 


1.171 1.168 1.150 1.117 



Deformation of the components (WINK model): 

a b c 4th ord. 5th ord. ro 

Star p 0.199 0.196 0.193 0.002 0.001 0.196 

Star s 0.221 0.217 0.212 0.004 0.001 0.217 





Primary 


Secondary 


Absolute dimensions: 






M/Mq 


11.28 ±0.46 


11.61 ±0.47 


R/Rq 


7.15 ±0.12 


7.92 ±0.13 


log g (cgs) 


3.782 ±0.023 


3.706 ±0.022 


Photometric data: 






V 


15.490 


15.318 


(b-y) 


-0.044 


-0.041 


mi 


0.062 


0.077 


a 


0.012 


0.027 


E( b -y) 


0.068 


0.082 


A v 


0.293 


0.309 


(b-y)o 


-0.112 


-0.111 


CO 


-0.002 


0.013 


[u-b] 


0.118 


0.054 


logT cff 


4.384 ± 0.005 


4.373 ± 0.005 


M bo i 


-5.75 ±0.06 


-5.86 ±0.07 


log L/Lq 


4.19 ±0.02 


4.23 ±0.03 


B.C. 


-2.31 


-2.25 


M v 


-3.44 ±0.06 


-3.61 ±0.07 


V - M v 


18.63 : 


± 0.08 


Distance (kpc) 


52.6 : 


± 2.0 



The combined intrinsic uvby indices given in Table 0] 
and the empirical [u — b] calibration by Napiwotzki et al. 
(1992| yield a mean T c s of 22000 K, and since the com- 
ponents have nearly identical surface fluxes, this value 
has been assigned to the primary component (eclipsed 
at phase 0.00) and kept fixed throughout the analysis. 
Linear limb darkening coefficients x (see Table I13|) by 
Diaz-Cordoves et al. (19HJ were adopted, and gravity 
darkening exponents of 0.25 (WINK) and 1.00 (WD), and 
bolometric albedos of 1.0, were assumed in accordance 
with atmospheres in radiative equilibrium. For the model 
calculation of the stellar deformations, a mass ratio of 
q = 1.0 was assumed, and the rotation of the stars was 
assumed to be pseudosynchronized (rate at periastron). 



7.3.1. WINK analyses 

Normally, at least i, ecosu, esinaj, r p , T s , and possibly k, 
would be kept free and determined through the differen- 
tial least squares procedure, but initial tests showed that 
this approach is not possible for HV12578, primarily be- 
cause i is very close to 90°. We have therefore calculated 
grids of solutions for adopted values of i, each of them 
for a range of adopted k values centered at 1.0. Fig. El 
illustrates the results from such WINK solutions for the y 
light curve. It is seen that the best fits are obtained for i 
between 89?8 and 90?0 and k between 0.98 and 1.02. For a 
given fc, most other parameters - especially the individual 
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Fig. 12. Results from WINK analyses of the y light curve of HV12578 for adopted i values between 89?7 and 90?0, 
and k = r s /r p values between 0.97 and 1.03. A constant phase shift of 0.00065 has been used. The individual figures 
show s: r.m.s. error (mag), e: orbital eccentricity e, w: longitude of periastron u) (°), rl: relative radius for primary 
component r p , r2: relative radius for secondary component r s , rl2: sum of relative radii, 121: luminosity ratio L s /L p , 
jn: surface flux ratio, t2: effective temperature of secondary component T s (K). Symbols are: cross i = 90?0, diamond 
i = 89?9, triangle i = 89?8, square i = 89?7. 



radii - are very little dependent on i within the range men- 
tioned above. Similar analyses for the other three bands 
support these findings, but are less accurate due to the 
fewer observations. 

Since i = 90? formally gives the best fit we adopt this 
value, and in Fig. ^| solutions from all four light curves 
are shown. The y and b results agree well, but due to the 
incomplete phase coverage in v and u, the elements from 
these two bands are less reliable. For the final solutions, 
we have chosen for all bands to fix e and uj at the y results. 
Such solution for k = 0.98, 1.00, and 1.02 are presented 
in Table ITTA and the O — C residuals for the k — 1.00 
solutions are shown in Fig. 1141 The uvby solutions agree 
quite well, and no systematic trends exist in the residuals. 
In all four bands, the component eclipsed near phase 0.49 
and by us referred to as the secondary, is slightly more 
luminous than the (primary) component eclipsed at phase 
0.00. 



7.3.2. Wilson-Devinney analyses 

For comparison, a few WD solutions for the y light curve, 
done for very nearly the same ratio of radii k as the 
WINK solutions listed in Table IT31 are given in Table [T4l 
Taking into account the effect from the different treatment 
of the expansion effects on the resulting radii, explained 
in Sect. Ifi~4"2l we find only small differences between the 
WINK and WD results. We notice a slightly higher surface 
flux ratio (via T s ) for the WD solution. 

7.3.3. Adopted photometric elements 

The adopted mean photometric elements are listed in 
TablelTSI The components are very nearly identical, and as 
seen, their relative sizes are slightly larger than for HV982. 
Relative radii accurate to 1% have been established, to- 
gether with well established elements for the very slightly 
eccentric orbit. Thus, when the results from the analy- 
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Fig. 13. Results from WINK analyses of HV12578 for an adopted orbital inclination i = 90?0, and k — r s /r p 
values between 0.95 and 1.05. The individual figures show s: r.m.s. error (mag), e: orbital eccentricity e, w: longitude 
of periastron ui (°), rl: relative radius for primary component r p , r2: relative radius for secondary component r s , 
rl2: sum of relative radii, 121: luminosity ratio L s /L p , jn: surface flux ratio, t2: effective temperature of secondary 
component T s (K). Symbols are: cross y, diamond 6, triangle v, square u. 



sis of the VLT/UVES spectra are available, we expect to 
be able to present accurate astrophysical parameters and 
distance for HV12578. 

8. HV1433 (SMC) 

At the time of the uvby observations, only about 30 SMC 
eclipsing binaries were known, and the selection of can- 
didates therefore more difficult than for LMC. Also, the 
quality of the photographic light curves from the Harvard 
survey (Gaposhkin I1977[l is rather low, and it was not 
possible to identify any clearly well-detached systems. 

CCD photometry has been obtained for several 
Harvard SMC systems (HV1620, Pritchard et al. I19i)8al 
HV1761, Duncan et al. 115531 HV1876, Jensen et al. 115551 
HV2016, Jensen et al. 115581 HV2208, W est et al. IT5521 
HV2226, Jensen et al. 115551 Bell et al. [T55TJ), and we 
initially included four additional in our project, HV1346, 
HV1433, HV1597, and HV11284. For HV1346 (P = 66 d 9) 



and HV1597 (P = 4.2), rather few observations were ob- 
tained, partial light curves were obtained for HV1433, 
and for HV11284 we obtained complete light curves. In 
the present paper, we will only include the HV1433 and 
HV11284 observations. 

As mentioned in Sect. [21 HV1433 is located in the 
southern part of the central region of SMC, just south of 
the OGLE - II SMC.SC4 field. Except for a 1.8 mag (y) 
fainter star about 3'/7 east of HV1433 (see Fig. OJ, which 
causes no problem for the PSF photometry, the surround- 
ing field is quite clean. 

Apart from the photographic Harvard light curve, no 
earlier observations of HV1433 have been published. 

8.1. Light curves 

HV1433 was observed on 19 nights between JD2448971 
and JD2450040, and 84(y), 38 (&), 33 (u), and 30 (u) points 
were obtained in the four bands. The light curves are given 
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Fig. 14. O - C residuals for the WINK solutions of HV12578 for i = 90° and k = 1.00 given in Table IT31 plotted 
versus orbital phase, seeing, sky background, and DoPHOT error (mag). 



in Tables 33-36 and shown in Fig. ^| Typical accuracies 
per point are given in Tabled 

Unfortunately, the coverage of the light curves is far 
from complete, and it is problematic to evaluate from the 
light curve shapes if HV1433 could maybe be detached, or 
if it is a semi-detached system, possibly with the secondary 
component filling its Roche lobe, like e.g. HV2208 and 
HV2226. However, the latter appears to be the most likely 
case. As the depth of secondary eclipse is only about half 
of that of primary, the temperatures of the components 
must be quite different. The orbit appears to be circular. 
In Sect. 18. 31 we will sketch possible photometric solutions. 



8.2. Ephemeris 

One time of primary minimum and two of secondary min- 
imum, derived from the uvby CCD observations, are given 
in Table IT5I together with the epoch of the ephemeris by 
Payne-Gaposchkin & Gaposchkin ( 19621) • We adopt the 



following ephemeris, derived from the new time of primary 
minimum and the old epoch 

Min I = 2450037.630 + 2^047002 x E , . 

±5 ±5 ' 

The period is slightly larger than that derived from the 
photographic Harvard observations (2^046987) . 

8.3. Light curve analysis, photometric elements 

We have applied only the Wilson-Devinney model for the 
photometric analysis of HV1433, since its components are 
expected to be larger and more deformed than appropri- 
ate for the WINK model; see Sect. 16.41 for details on the 
actual WD code. A semi-detached configuration with the 
secondary filling its Roche lobe has been assumed, since 
initial tests showed that detached models converges to- 
wards this configuration. 

The standard uvby indices for the combined light of the 
components, listed in Table 01 correspond to an effective 
temperature of about 23000 K (Napiwotzki et al. Q[992), 
which has then been adopted for the primary component, 
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Fig. 15. uvby light curves of HV1433 (SMC). The theoretical curves correspond to the solutions (q = 0.50) given in 
Table IT71 



Table 12. Times of minima for HV12578 (LMC). 
Phases are calculated from the ephemeris given in Eq. [3] 
References are: G1977, Ganoshkin IWF1 C2003, this pa- 
per. The r.m.s. errors of the G1977 times are estimated. 
The minimum at JD2424501.614 listed by G1977 has been 
omitted. 



HJD-2400000 


r.m.s. 


Type 


Phase 


Reference 


16823.638 


0.050 


P 


0.010 


G1977 


23681.879 


0.050 


P 


0.010 




23907.499 


0.050 


P 


-0.028 




26303.577 


0.050 


P 


-0.015 




26328.446 


0.050 


P 


0.012 




27800.507 


0.050 


s 


0.494 




29881.575 


0.050 


s 


0.507 




29927.447 


0.050 


p 


0.001 




29994.482 


0.050 


p 


0.027 




30640.598 


0.050 


s 


0.518 




31496.277 


0.050 


s 


0.498 




48971.777 


0.001 


p 


-0.000 


C2003 


49343.835 


0.001 


p 


0.000 




49637.738 


0.001 


s 


0.492 




50040.825 


0.010 


p 


0.002 





Table 14. WD solutions (y) for HV12578 for fixed Sl p . 
i = 90?0, T p = 22000if and x = 0.32 have been assumed. 





y 


y 


y 


e 


0.0157 


0.0153 


0.0171 




±12 


±14 


±20 


o 


150.7 


150.3 


139.4 


±6.9 


±8.0 


±6.9 




5.040 


5.120 


5.200 


n s 


5.117 


5.126 


5.138 




±18 


±16 


±17 


r p 


0.2513 


0.2462 


0.2414 


r a 


0.2464 


0.2458 


0.2452 


k 


0.981 


0.998 


1.016 


r p + r a 


0.4977 


0.4920 


0.4866 


Ts (K) 


22502 


22512 


22461 


±54 


±57 


±57 


L s f L p 


0.983 


1.020 


1.052 


cr(mag) 


0.0115 


0.0114 


0.0121 



although it might be slightly on the low side due to the 
effect on the indices of the cooler component. Linear limb 
darkening coefficients by Diaz-Cordoves et al. I)1995|l were 
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Table 17. WD solutions for HV1433 for adopted mass ratios between 0.4 and 0.6. A semi-detached configuration is 
assumed, and T p — 23000 K has been adopted throughout. 





y 


y 


L 




V 


u 


y 


q — m s /m p 


0.40 


0.50 


0.50 


0.50 


0.50 


0.60 




88.37 


86.17 


84.43 


86.15 


86.68 


84.80 


±20 


±15 


±18 


±37 


±44 


±10 


O p 


3.3394 


3.4093 


3.5562 


3.4710 


3.4474 


3.6630 




±84 


±192 


±398 


±557 


±425 


±201 


Q, 


2.6781 


2.8758 


2.8758 


2.8758 


2.8758 


3.0624 


r p 


0.347 


0.352 


0.334 


0.344 


0.347 


0.333 


r s 


0.303 


0.321 


0.321 


0.321 


0.321 


0.336 


X p 


0.34 


0.34 


0.38 


0.40 


0.40 


0.34 


Xs 


0.37 


0.37 


0.43 


0.46 


0.43 


0.37 


T s (K) 


16157 


16150 


16425 


16577 


16595 


16180 


±67 


±54 


±70 


±102 


±117 


±51 


L s 1 ' Lp 


0.455 


0.494 


0.539 


0.494 


0.356 


0.600 


a(mag) 


0.0134 


0.0111 


0.0098 


0.0148 


0.0150 


0.0120 



Table 15. Adopted photometric elements for HV12578. 
r p and r s are corrected for expansion due to rotation and 
tidal effects, the r$ result (both components) is not. The 
individual flux and luminosity ratios are based on the 
mean stellar and orbital parameters. 



i 


90?0±0?2 


ecosu; 


-0.01341 ± 0.00030 


esinw 


0.00745 ± 0.00200 


e 


0.0153 ±0.0010 


UJ 


151?0±6?5 


r p 


0.2425 ± 0.0025 


r s 


0.2425 ± 0.0025 


r p + r 3 


0.4850 ± 0.0035 


k 


1.00 ± 0.02 



y b v u 

Js/Jp 1.027 1.026 1.039 1.011 
L s /L p 1.028 1.027 1.040 1.013 

Deformation of the components (both; WINK model): 

a b c 4th ord. 5th ord. ro 

0.2498 0.2442 0.2408 0.0037 0.0009 0.2449 



Table 16. Times of minima for HV1433 (LMC). 
Phases are calculated from the ephemeris given in Eq. 
References are: SA09, Payne-Gaposchkin & Gaposchkin 
{HSU; C2003, this paper. 



HJD-2400000 


r.m.s. 


Type 


Phase 


Reference 


24431.682 




P 


0.0000 


SA09 


50037.6300 


0.0050 


P 


0.0000 


C2003 


50038.6544 


0.0005 


S 


0.5004 




50040.7055 


0.0010 


s 


0.5024 





8.4. Discussion 

Due to its brightness, its position within the SMC, and 
its uncomplicated surrounding field, HV1433 is an inter- 
esting target for future photometric and spectroscopic ob- 
servations and studies, primarily for distance determina- 
tions. On the other hand, many of the new eclipsing bi- 
naries in the central region of SMC, now known from the 
OGLE — II project, could just as well be considered. They 
include about 150 detached systems, most of which are, 
however, much fainter than HV1433 and therefore in gen- 
eral more difficult to observe accurately enough. 



adopted together with gravity exponents and reflection 
albedos of 1.0. A circular orbit and synchronized rotation 
of the components were assumed. 

Solutions for a range of fixed mass ratios indicate that 
values below 0.4 and above 0.6 are not very likely. In 
Table 1171 we present a few solutions within the remain- 
ing range for which acceptable fits are reached. Definitive 
solutions can not be obtained from the incomplete light 
curves, but relative radii of around 0.34 (primary) and 
0.32 (secondary), a temperature difference between the 
components of about 7000 K, and an orbital inclination 
of about 86° arc estimated. 



9. HV11284 (SMC) 

As mentioned in Sect. [21 HV11284 is located slightly 
west of the optical center of SMC in the OGLE — II field 
SMC.SC5 and is identical to the OGLE - II eclipsing bi- 
nary SMC_SC5_140701. We selected HV11284 as one of 
the most promising systems among the few Harvard can- 
didates available and obtained complete CCD light curves 
1992-94. In addition, BVI light curves from 199 7-98 ar e 
now also available from OGLE — II (Udalski et al.E?98b). 

Being nearer to the SMC center, the HV11284 field is 
more crowded than that of HV1433 with a higher density 
of faint stars. There are no bright, close companions to 
HV11284, but the shape of its image on CCD exposures 
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Fig. 16. uvby light curves of HV11284 (SMC). The theoretical curves correspond to the solutions given in Table l2"ll 



obtained at about 1 arcsec seeing indicates the presence 
of a possible very close, but also quite faint companion 
slightly to the east; see Fig.0] If real, it is included in the 
DoPHOT photometry of HV11284 and therefore may add 
a small amount of 3rd light to the light curve data. We 
estimate that it is at least 3-4 mag fainter than HV11284 
at maximum light, leading to an upper level of about 6% 
of possible 3rd light (L 3 /(L p + L s )). 



9.1. Light curves 

HV11284 was observed on 34 nights between JD2448971 
and JD2449638, and 219 (y), 81 (6), 70 («), and 63 (u) 
points were obtained in the four bands. The light curves 
are given m Tables 37-40 and shown in Fig. [TBI Typical 
accuracies per point are given in Table El 

As seen from the relatively broad eclipses and the light 
variations outside eclipses, HV11284 consists of two rather 
close and significantly deformed components. It is, how- 
ever, not quite obvious if the system is detached or semi- 
detached. The eclipses differ significantly in depths, and 
the secondary component eclipsed at phase 0.5 is clearly 



Table 18. Times of minima for HV11284 (LMC). 
Phases are calculated from the ephemeris given in Eq. |3J 
References are: SA09, Payne-Gaposchkin & Gaposchkin 
(UJSJ; U1998, Udalski et al. I|1998b|> adjusted by 
-0<?0087; C2003, this paper. 



HJD-2400000 


r.m.s. 


Type 


Phase 


Reference 


34682.321 




P 


0.0000 


SA09 


48975.607 


0.005 


S 


0.4971 


C2003 


49343.598 


0.002 


P 


-0.0004 




49633.636 


0.002 


P 


0.0003 




50471.1116 




P 


0.0001 


U1998 



much cooler than the primary. The orbit appears to be 
circular. 



9.2. Ephemeris 

Two times of primary minimum and one of secondary min- 
imum, derived from the uvby CCD observations, are given 
in Table UhI together with the epochs of the ephemeris by 
Payne-Gaposchkin & Gaposchkin ( 196SJ) and Udalski et 
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Table 19. WD solutions for HV11284 based on the two 
available OGLE / light curves; see text for details. A mass 
ratio q = 0.80 and T p = 23000 K has been adopted. 





I(DIA) 


/ 


if) 


84.32 


84.66 


±14 


±44 




3.6172 


3.6325 




±49 


±138 


o s 


4.1145 


4.0562 




±97 


±372 


r p 


0.367 


0.365 


r s 


0.270 


0.275 


X p 


0.24 


0.24 


Xs 


0.31 


0.31 


Ts (K) 


15242 


15183 


±49 


±95 


L s 1 Lp 


0.295 


0.309 


a (mag) 


0.0090 


0.0171 



al. H1998b(l . Based on a reanalysis of the OGLE / light 
curve, the latter has been adjusted by — O^OOS?. We note 
that the light curve tabulated by Gaposhkin 1)1977(1 is 
shifted with primary minimum located at phase 0.25, i.e. 
as if the epoch should correspond to maximum light level. 
From a weighted least squares fit to the times of primary 
eclipse and the two epoch, we derive the following linear 
ephcmeris: 

Min I = 2449343.5996 + 3 d 6254398 x E frS 
±8 ±4 (5) 

The period is somewhat shorter than that determined 
from the photographic Harvard observations (3^625487), 
and somewhat longer than the result (3*?62537) published 
by Udalski et al. ; 1 99M > i . However, a reanalysis of their I 
(DIA; see below) light curve yields 3^62545. 

9.3. Light curve analysis, photometric elements 

From the shape of the light curves, the two components 
of HV11284 are expected to be larger and more deformed 
than appropriate for the WINK model, and we have there- 
fore applied only the Wilson-Devinney model for the pho- 
tometric analysis; see Sect. 16.41 for details on the actual 
WD code. The standard uvby indices for the combined 
light of the components, listed in Tabled! correspond to an 
effective temperature of about 23000 K (Napiwotzki et al. 
1992), which has then been adopted for the primary com- 
ponent. As for HV1433, the actual effective temperature 
might be slightly higher, but tests show that the possible 
difference is of no importance for the light curve analysis. 
Linear limb darkening coefficients by Diaz-Cordoves et al. 
(1995 ) were adopted together with gravity exponents and 
reflection albedos of 1.0. A circular orbit and synchronized 
rotation of the components were assumed. 

The stellar shapes depend on the mass ratio between 
the components, and on their rotational velocities, which 



Table 20. WD solutions for HV11284 (y light curve) 
for adopted mass ratios between 0.6 and 1.4. A detached 
configuration is assumed, and T p = 23000 K has been 
adopted throughout. 



q = m 3 /nip 


0.60 


0.80 


1.00 


1.20 


1.40 


iO 


84.78 


84.95 


85.37 


84.38 


82.75 


±11 


±10 


±8 


±10 


±7 


Op 


3.3971 


3.6418 


3.9088 


4.1889 


4.4898 




±61 


±49 


±30 


±39 


±42 


a, 


3.3725 


4.0714 


4.8617 


5.5874 


5.9983 




±59 


±101 


±94 


±148 


±175 


r p 


0.3679 


0.3633 


0.3557 


0.3466 


0.3353 


r a 


0.2824 


0.2737 


0.2622 


0.2588 


0.2726 


F P 


0.91 


0.94 


0.91 


0.97 


0.95 


F s 


0.84 


0.82 


0.74 


0.72 


0.72 


dp 


0.34 


0.34 


0.34 


0.34 


0.34 


Xs 


0.37 


0.37 


0.37 


0.37 


0.37 


Ts (K) 


15160 


14790 


14393 


14212 


14249 


±40 


±37 


±35 


±34 


±34 


L s 1 ' Lp 


0.314 


0.286 


0.260 


0.259 


0.304 


a (mag) 


0.0097 


0.0086 


0.0080 


0.0078 


0.0079 



Table 21. WD solutions for HV11284 (OGLE I (DIA) 
light curve) for adopted mass ratios between 0.6 and 1.4. 
A detached configuration is assumed, and T p = 23000 K 
has been adopted throughout. 



q = m s /m p 


0.60 


0.80 


1.00 


1.20 


1.40 


in 


83.22 


84.32 


84.18 


83.29 


81.86 


±23 


±14 


±14 


±9 


±7 




3.3861 


3.6172 


3.8983 


4.1957 


4.5005 




±97 


±51 


±47 


±35 


±41 


n s 


3.3861 


4.1114 


4.8585 


5.4805 


5.8599 




±142 


±97 


±184 


±128 


±122 


r p 


0.3695 


0.3668 


0.3572 


0.3457 


0.3340 


r a 


0.2868 


0.2696 


0.2625 


0.2651 


0.2803 


J p 


0.24 


0.24 


0.24 


0.24 


0.24 


X s 


0.31 


0.31 


0.31 


0.31 


0.31 


Ts (K) 


15648 


15242 


14999 


14870 


14985 


±58 


±49 


±45 


±43 


±43 


L s / L p 


0.349 


0.294 


0.279 


0.302 


0.364 


cr(mag) 


0.0100 


0.0090 


0.0080 


0.0075 


0.0075 



are unfortunately not known since spectroscopic obser- 
vations are not available. This means that we needed to 
explore various plausible configurations. 

Besides the uvby light curves, we have also analyzed 
the / light curve available from the OGLE WWW cata- 
logue 2 . B and V light curves are also available, but they 
contain too few points for meaningful analyses. Two ver- 
sions of the light curves are available: The original light 
curves are based on photometry derived from a modi- 
fied DoPHOT package, and recently new and more accu- 
rate light curves based on difference image analyses (DIA) 
have been published (Zebrun, Soszynski, & Wozniak 2001 
Zebrun et al l2"tj0l")l . We have used the I (DIA) light curve 
for the WD analyses. It contains 312 points obtained 

2 http://sirius.astrouw.edu.pl 
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Table 22. WD solutions for HV11284 (y light curve) 
for adopted mass ratios between 0.6 and 1.4. A detached 
configuration has been assumed, and T p = 23000 K has 
been adopted throughout. A small amount of 3rd light of 
0.06 (L%/(Lp + L s )) has been included. 



q = m s /m p 


0.60 


0.80 


1.00 


1.20 


1.40 


in 


86.02 


86.04 


86.06 


86.29 


84.83 


±13 


±13 


±10 


±8 


±8 




3.4103 


3.6442 


3.9034 


4.1715 


4.4576 




±47 


±45 


±35 


±34 


±40 


n s 


3.2994 


3.9514 


4.6699 


5.4515 


6.0543 




±42 


±83 


±80 


±96 


±93 


r p 


0.3659 


0.3629 


0.3565 


0.3490 


0.3394 


r s 


0.2931 


0.2859 


0.2766 


0.2669 


0.2696 


X p 


0.34 


0.34 


0.34 


0.34 


0.34 


X s 


0.37 


0.37 


0.37 


0.37 


0.37 


T s 


15280 


14934 


14595 


14307 


14175 




±38 


±35 


±33 


±33 


±34 


L s 1 Lp 


0.347 


0.318 


0.293 


0.273 


0.288 


a (mag) 


0.0095 


0.0085 


0.0077 


0.0075 


0.0079 



Table 24. WD solutions for HV11284 for an adopted 
mass ratios of q = 0.4. A semidetached configuration with 
the less massive secondary component filling its Roche 
lobe has been assumed, and T p = 23000 K has been 
adopted throughout. 





I(DIA) 


y 


b 


V 


u 




81.77 


82.71 


82.26 


82.10 


82.97 


±6 


±5 


±8 


±8 


±16 




3.4145 


3.3410 


3.4436 


3.7161 


3.3378 




±142 


±81 


±110 


±105 


±253 


n s 


2.6781 


2.6781 


2.6781 


2.6781 


2.6781 


r p 


0.338 


0.347 


0.334 


0.324 


0.347 


r a 


0.303 


0.303 


0.303 


0.303 


0.303 


X p 


0.24 


0.34 


0.38 


0.40 


0.40 


X s 


0.31 


0.37 


0.43 


0.46 


0.43 


Ts (K) 


16056 


15638 


15773 


16175 


15166 


±43 


±37 


±85 


±61 


±70 


L s j L p 


0.506 


0.446 


0.462 


0.486 


0.305 


cr(mag) 


0.0081 


0.0092 


0.0085 


0.0086 


0.0119 



Table 23. WD solutions for HV11284 for an adopted 
mass ratio of q = 0.8. A detached configuration has been 
assumed, and T p — 23000 K has been adopted throughout. 





y 


b 


V 


u 




84.95 


85.00 


85.10 


84.52 


±10 


±20 


±32 


±45 


O p 


3.6418 


3.6849 


3.7161 


3.7443 




±49 


±76 


±105 


±162 


O s 


4.0714 


4.1212 


4.1270 


3.9734 




±101 


±164 


±288 


±479 


r P 


0.363 


0.357 


0.353 


0.349 


r s 


0.274 


0.269 


0.268 


0.284 


X p 


0.34 


0.38 


0.40 


0.40 




0.37 


0.43 


0.46 


0.43 


Ts (K) 


14790 


14852 


15310 


14178 


±37 


±65 


±72 


±79 


L s 1 Lp 


0.286 


0.272 


0.275 


0.212 


a (mag) 


0.0086 


0.0086 


0.0084 


0.0101 



through a large number of nights. It is generally well co- 
vered, but rather few points are obtained during primary 
eclipse. The original / light curve leads to nearly identical 
but less well defined WD solutions; see Table 

First, we have considered HV11284 to consist of two 
detached components, and a set of y and / solutions for 
a range of adopted mass ratios between 0.6 and 1.4 are 
listed in Table I2UI and 1211 For all mass ratios, the hotter 
primary component is larger than the 8-9000 K cooler se- 
condary component. Filling factors F (see e.g. Mochnacki 
1984J are between 0.91 and 0.97 for the primary, and be- 
tween 0.84 and 0.72 for the secondary. Formally, the best 
solutions correspond to a mass ratio of around 1.2, which 
is, however, not very likely. For mass ratios above 1.0, and 
given the large temperature difference, the more massive 
component is expected to have evolved off the main se- 
quence and to have become not only the cooler but also 



the larger component. This is clearly not supported by 
the light curve solutions. If, on the other hand, both com- 
ponents of HV11284 are still within the main sequence 
band, current evolutionary models predict a mass ratio 
not larger than about 0.6 for a primary temperature of 
around 23000 K and a 8-9000 K temperature difference 
between the components. Attempts to improve the fit of 
the theoretical light curves to the observations for mass 
ratios around 0.6 by adjusting rotational velocities, reflec- 
tion coefficients, and limb darkening coefficients (includ- 
ing using non-linear laws) have not been very successful, 
and adding a small amount of third light from the pos- 
sible close companion does not help much either, as seen 
from Table|221 Without knowing the mass ratio, it is, how- 
ever, still considered too early to exclude definitively that 
HV11284 is a detached system. Analyses of the b, v, and 
u light curves lead to similar but less accurate results and 
conclusions. In Table l2*3l we compare the uvby solutions 
for an arbitrary (although not very likely) mass ratio of 
0.8 which, at the moment, seem to represent the deforma- 
tion of the components adequately, although not as well 
as higher mass ratios. 

Next, we have considered HV11284 to be a semi- 
detached system with the cooler component filling its 
Roche lobe. For this configuration, we also find it impossi- 
ble to determine safely the mass ratio from the light curve 
solutions. It is, however, seen that only values below 0.5- 
0.6 lead to acceptable agreement between the theoretical 
light curves and the observations. The quality of the fits is 
comparable to those for detached solutions, meaning that 
a semi-detached configuration for HV11284 is certainly 
possible. Solutions for an adopted mass ratio of q = 0.4 
are listed in Table |2H 

In conclusion, we find relative radii of 0.33-0.36 (pri- 
mary) and 0.26-0.30 (secondary) for the components and 
a temperature difference of of 8-9000 K between them, 
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but the actual configuration of HV11284 remains an open 
question. The orbital inclination is 82°-85°. 

9.4. Discussion 

The uvby and OGLE I light curves of HV11284 are well 
covered and of quite high quality. They should allow defini- 
tive photometric solutions once the mass ratio is estab- 
lished from spectroscopic observations. HV11284 is poten- 
tially a valuable system for distance determination, and 
since it is relatively bright, we suggest that such observa- 
tions are done, e.g. at the VLT. 

10. Concluding remarks 

With present-day facilities, precise studies of eclipsing bi- 
naries in the LMC and the SMC have become possible. 
Such objects are excellent distance indicators, and dis- 
tance moduli better than ±0.10 mag can often be ob- 
tained. Also, accurate absolute dimensions for MC binary 
components are important for the study of internal struc- 
ture, mass loss and evolution of massive metal-deficient 
stars. 

We have in this paper, as part of a large scale project 
on MC eclipsing binaries, presented new complete uvby 
light curves and analyses of them for the LMC systems 
HV982 and HV12578, and the SMC system HV11284, 
and partial results for the SMC system HV1433. Accurate 
standard uvby indices have also been established for each 
binary, and individual interstellar reddenings have been 
determined. 

HV982 and HV12578 are well-detached systems with 
eccentric orbits, each consisting of components of compa- 
rable sizes, and accurate photometric elements have been 
established. Adopting the spectroscopic elements given by 
Fitzpatrick et al. (2002J for HV982, we have derived ab- 
solute dimensions of its component which agree well with 
their results. A distance modulus of V$ — My = 18.63±0.08 
is obtained, corresponding to a distance of 52.6 ± 2.0 
kpc, which is in formal agreement with (although slightly 
larger than) their determination. In a forthcoming paper, 
we will combine our photometric results with those un- 
derway from high-dispersion spectra obtained with the 
UVES spectrograph on the ESO VLT and present abso- 
lute dimensions, chemical abundances and distances for 
HV12578 and refined results for HV982. 

HV1433 and HV11284 (identical to the OGLE - II sys- 
tem SMC_SC5_140701) both consist of two rather close, 
deformed and quite different stars. Because the mass ra- 
tios between the components - and their rotation rates 
- are not known, definitive photometric elements could 
not be obtained. We have, however, presented a sample 
of possible photometric solutions. Since HV11284 is re- 
latively bright and has complete uvby I light curves, we 
recommend that spectroscopy is obtained for this SMC 
system. In addition to spectroscopy, more photometry is 
also needed for HV1433. 
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